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We investigate the isotropic two-leg S = \ ladder with gen- 
eral bilinear and biquadratic exchange interactions between 
spins on neighboring rungs, and determine the Hamiltonians 
which have a matrix product wavefunction as exact ground 
state. We demonstrate that a smooth change of parameters 
leads one from the S = i dimer and Majumdar-Ghosh chains 
to the S = 1 chain with biquadratic exchange. This proves 
that these model systems are in the same phase. 
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Low dimensional quantum antiferromagnets have at- 
tracted a large amount of both theoretical and experi- 
mental interest in recent years. Theoretically, the results 
of calculations for model systems in one dimension (ID) 
clearly show the difference between systems with a gap- 
less spectrum of excitations and power law decay of spin 
correlations and gapped systems with exponentially de- 
caying correlation functions. The main representatives of 
these two classes are resp. the S = \ chain with isotropic 
nearest neighbor (NN) exchange interaction jl| and the 
isotropic £ = 1 (Haldane) chain Q . Gapped elementary 
excitation spectra are also found in more complicated 
low-dimensional S — g systems such as chains with suf- 
ficiently strong additional next-nearest neighbor (NNN) 
exchange or with alternating exchange (with noninteract- 
ing dimers as the simplest limit) and also isotropic spin 
ladders with an even number of legs (^] . 

Many of these models are realized to a high degree of 
accuracy in simple compounds as demonstrated by the 
following examples: KCUF3 (ID isotropic S = \ anti- 
ferromagnet [||), CaCuGe 2 Og (weakly interacting dimers 
@), (V0 2 )P 2 7 |] and SrCu 2 3 @ (two leg spin ladder) 
and Ni(C 2 H 8 N2)2N0 2 C104 (NENP, ID S = 1 antiferro- 
magnet ||). Although the precise spin Hamiltonians of 
these substances are not known, the study of the simple 
theoretical models mentioned above has contributed in 
an essential way to the understanding of the behaviour 
of the real materials and these are believed to be in the 
same phase as the simple models. 

In recent years arguments and numerical evidence have 
been presented that the various models with a gapped ex- 



citation spectrum, in particular the S = | dimer chain, 
the gapped nearest and next nearest neighbor exchange 
chain, the two-legged S = \ spin ladder and the S = 1 
(Haldane) chain are all in the same phase §-@. In this 
letter we present the following proof for this conjecture: 
We describe a class of Hamiltonians with exactly known 
nondegenerate ground states and show that it is possi- 
ble, by a continuous change of parameters within this 
class, to connect smoothly the hamiltonians of the ID 
dimer chain, the Majumdar-Ghosh chain, a generalized 
spin ladder (which includes biquadratic interactions) and 
of the S = 1 antiferromagnetic chain with additional bi- 
quadratic exchange. Given the relevance of these model 
systems with exactly known ground states for the tra- 
ditional model systems and for the real quasi ID com- 
pounds this proves the existence of a single phase for 
the gapped low dimensional spin systems above. As a 
"byproduct" of this investigation, we also present a new 
set of models of frustrated S — \ spin chains (including 
only bilinear NN and NNN interactions) whose ground 
states can be found exactly. 

Our approach starts from the observation [jl3| that the 
exactly known ground states of the Shastry-Suthcrland 
S= \ chain ^| 

H = ]T(1 + (-l)"5)S n S n+1 + \{l - 6)S n S n+2 , (1) 

n 

which includes the Majumdar-Ghosh (S = 0) jl6| and 
dimer (S — 1) limits as special cases, and the 5=1 
chain with biquadratic exchange 

H = S„S„+i — /3(S„S„+i) 2 (2) 

n 

at /? = -i (the AKLT model which is believed to 
be in the same phase as the Haldane chain with ft = 0) 
have very similar structure: both ground states can be 
written as matrix product (MP) wave functions, 

I*.) -lift, ft-^V5a|t_0i b\s) t -a\t Q )J- (3) 

The ground state of the Shastry-Sutherland chain is ob- 
tained setting a — b — i, where |s)j and 1 1)^ are the sin- 
glet, resp. triplet formed by two spins S = 5 on adjacent 
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sites n and n + 1 (n even for 5 > in (0)). The ground 
state of the AKLT chain is obtained by setting a = -j- , 
6 = 0, and identifying the remaining triplet with the three 
states of the spin S — 1 at site i. The matrix elements 
of the MP wavefunction l^o) are related to boundary ef- 
fects for open boundary conditions whereas for periodic 
boundary conditions the trace should be taken jl7| . Ma- 
trix product (or finitely correlated) wave functions were 
first introduced by Fannes et al [fL8l a decade ago and 
have since then found widespread applications in exact 
and variational calculations |l9|-[2^]. These wave func- 
tions are characterized by finite (and typically rather 
short) correlation lengths and are therefore tailored to 
deal with gapped systems. 

It should be emphasized that our MP ansatz for the 
ladder (||), in contrast to that used in Ref. [p2| , respects 
rotational symmetry of the Hamiltonian (M), i.e. the 
wavefunction \tp ) under periodic boundary conditions is 
a global singlet for any value of the parameter u |l3],^4| . 
The ansatz used by Su has a 'built-in' broken rota- 
tional symmetry and thus can access only states of fer- 
romagnetic type. 

We generalize the models of Eqs. (0,0) by considering 
the following spin ladder Hamiltonian (see Fig. [l]): 
H = J2i where 

h>i,i+i = Co + — ^oSi,iS2,i + — JoSi^+iSa.-j+i 

+JlSl 1 jSi ) j+l +^iS2,iS2,-t+l + ^2Sl,i S2,i+1 + ^^S.jSi^+i 
+-K"(Sl,iSi ) i+l)(S2,-iS2,i+l)+ K (81,182,1+1) (S2,iSi,i+i) 

+A" / (Si, i Sa,i)(Si, i+ iS a ,t + i) (4) 

This is the general isotropic S = ^ Hamiltonian with 
exchange interactions restricted to neighboring rungs of 
the ladder. For periodic boundary conditions 5 (Jo + J' a ) 
is the coupling on the rungs, in addition there are 4 pair 
exchange couplings on the legs and on the diagonals and 
3 biquadratic terms. Together with the constant Co we 
have a total of ten parameters. One combination of these 
is irrelevant since it sets the energy scale. It is essential 
to include the diagonal interactions as first introduced 
by White fL4j ) in order to obtain the various limits of 
interest. These limits include the isotropic spin ladder 
and also the S = 1 chain when strongly ferromagnetic 
exchange on the rungs (Jo — > —00) is considered to form 
5 = 1 units with an effective pair interaction J e ff and 
effective biquadratic interaction K e ff 

Jeff = \{Jl + J[ + J 2 + J 2 ) + \{K + K'), 



Keff = \{K + K') 



(5) 



Our procedure in the following is based on the ideas 
presented in |2(J and generalizes the work for the 5=1 
chain with biquadratic exchange |l7||: We start from the 




FIG. 0. A generalized spin-^ ladder with all possible isotropic 
exchange interactions between nearest-neighbor rungs. Ovals 
show spin pairs whose states are involved in the elementary 
matrix gi defined by (|^). Couplings are described in the text, 
see (§). 

MP wave function ([|), considering \s) and as the 
states of a single rung, and require that the parameters 
in the Hamiltonian (||) and the free parameter u = b/a 
in |i/>o) satisfy the following conditions: 
(i) |^o) is annihilated by H, i.e. 



gigi+i — ; 



(G) 



(ii) all other states of H have the energy E > 0. 
Then |-0o) is the ground state of H , with the energy den- 
sity —Co per rung. 

To proceed we write the local Hamiltonian h in the al- 
ternative formulation |21| using projectors on states with 
fixed angular momentum of the 4 spin plaquette + 1) 
of the ladder: 



£ = A 2 $>2,m><*2,m|+ £ A w Xo><*o,bl 

fej=l,2 



M 



fc,Z = l,2,3 



M 



(7) 



Here l^j^f ) are the multiplets which can be formed from 
the 4 spins S — ^ of one plaquette: one quintuplet, three 

(k I) (I k) 

triplets and two singlets; the 10 parameters Xj' ' = Ay 
are linearly related to the 10 parameters in the Hamilto- 
nian of Eq. (|). The multiplets |*^) are given in terms 
of the singlets |s) and triplets \t) on rungs i and i + the 
plaquette singlets are a linear combination of \ss) and 
\(tt)j = o), the plaquette triplets are linear combinations 
of \ts), \st) and \(tt)j = i) and the plaquette quintuplet is 
\{tt)j=2) ■ It is easily verified that only one singlet 



(3 + u 4 )- 1/2 (u 2 \ ss )-V3\tt) J=() ) 
and one triplet 



- I,T,(1)\ 



=0 )= $ 



(8) 



(1 + n 2 )- 1/2 {^(N + I st )) - I**>j=i} = ( 9 ) 

occur in gig^i- A convenient choice for the remaining 
multiplets is 
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|*£ 2) ) = (3 + u 4 )- 1 / 2 (y/3\ S3 ) + u 2 \tt) J=0 ) 
\^) = (1/V2)(|£ S ) - 1st)) 

\9<®) = (1 + u 2 )^ 2 {(l/V2)(\ts) + \st)) + u\tt) J=1 } 

1*2) = |tt),7=2 (10) 

Conditions (i) result in 5 equations, corresponding to 
A^ 1 ) = A^ 2) = A< ia) = Ai 1 ' 2 ' = Ai 1 ' 3 ' = 0. (11) 

Thus, in order to lead to MP ground states, the Hamilto- 
nian has to project on the space |* 2) ), |*i 3) ), 1*2) 

and condition (ii) results in the inequalities 

\ { 2) > 0, \{ a) > 0, A 2 > (12) 

where (a = 1, 2) are the two eigenvalues of the matrix 
•M > ihj = 2,3). The inequalities ( |l2] ) guarantee that 
the ground state is nondegenerate (apart from the fact 
that Eq. (||) describes four ground states with different 
behaviour at the boundaries of open ladders), since it can 
be shown by induction with respect to the ladder length 
that the plaquette states of Eqs. (||,||) do not allow a zero 
energy state which is different from the MP ground state. 

For an explicit discussion of the results we use a some- 
what simplified Hamiltonian with less freedom than in 
the general case: First, we set Jo = Jq, and from 
one of the conditions (i) we find that this necessarily 
means J± = J{ (which respects the ladder symmetry) and 
A^ 1,2 -* = 0; thus the eigenvalues A^ are identical to the 
Af ,2) and \f ' 3) . Further we set K-K' = K" = (these 
combinations are found irrelevant in a more detailed 
treatment which will be published separately). Then we 
are left with the 6 coefficients Co, Jo, Ji, J2, J'%, K and 
the conditions (i) take the following form 

2Ji - J 2 - J 2 + u 2 (4C /3 - J + K/2) = 

2(2 Ji + J 2 + J 2 ) - 4C - J - 7K/2 (13) 

W(J 2 + J' 2 - 2Ji) = 
2(J 2 - J 2 ) + u (4C - J + 2Ji -J2-J2- K/2) = 
2Ji + J 2 + J 2 + 3K/2 - 4C - Jo + 2u(J 2 - J' 2 ) = 

The general solution of these 4 linear equations contains 
2 arbitrary constants, when the parameter u = b/a of 
the MP wavefunction is fixed. We absorb one of the two 
constants in the energy scale, denote the remaining one 
by x and define 

u 2 = F 

to obtain the following family of hamiltonians with exact 
MP eigenstate: 

C* = 3 (9 + 2x - 3(1 + x)F + 3xF 2 ) /32 , 
Jo = 3 (2 - x - (1 + x)F + xF 2 ) /4 , 



2Ji + J 2 + J 2 = (15 - 9(1 + x)F + 7xF 2 ) /4 , 

2 Ji - J 2 - J 2 = -xF , 

2J 2 - 2J^ = F 1/2 (3(1 + x)/2- xF) , 

K = 3(F- l)(xF- l)/4. (14) 

With these definitions the conditions of Eq. (^2|) that the 
MP state is a ground state are obtained as follows 

A = x(3 + F 2 (x))/4 > 

A^ = (3(1 + x) + 2xF{x)) /8 > 

\ {2) = (1 + F(x)) (3(1 +x) - 2xF(x)/8) > 

A 2 = (18 + 8xF 2 (x) - 9(1 + x)F(x)) /8 > (15) 

The eigenvalues Ao and A^ are seen to be positive for 

— (2) 

x > 0, F > 0. The eigenvalues A^ ; and A 2 require a 
more detailed discussion; it is easily seen that a sufficient 
(though not necessary) condition for A 2 > is J e // > 0, 
i.e. an effective antiferromagnetic interaction between 
the two rungs of one plaquette, and the condition on A^ 
is equivalent to the requirement u( J 2 — J 2 ) > 0. It is con- 
venient to discuss the possible hamiltonians considering 
F as function of x. We use F = u 2 instead of u keeping 
in mind that changing the sign u — > — u amounts just to 
interchanging the ladder legs and thus does not bring in 
any new physics; therefore from now on we assume that 
F 1 ' 2 > 0. 

The following illustrative members of this family of 
hamiltonians are now easily obtained: 

(i) Shastry-Sutherland model: 

F=l, x = 36/(2 + 5), 0<<5<1. (16) 

Since the MP wave function with F = 1, i.e. a = b 
corresponds to singlets on one type of diagonals jl3) , this 
also covers (after a translation of one of the ladder legs) 
the case a = 0. The eigenvalue Ao vanishes for a; = 0, 
which is to be expected since in the Majumdar-Ghosh 
limit a second degenerate eigenstate exists (singlets on 
the alternative bonds). 

The solution (16) also applies for partially ferromag- 
netic interactions, x > 1, resp. (1 — S)/(l + S) < 0. 
Although the effective interaction between units on the 
diagonals is ferromagnetic, the Shastry-Sutherland dimcr 
state is the ground state up to x — 9 (corresponding to 
S = —3), since all As are positive. At the singular point 
5 = — 3 the eigenvalue A 2 vanishes and the energy of the 
dimer ground state coincides with the energy of the fully 
polarized ferromagnetic configuration. Thus this is the 
point of the first-order quantum phase transition from 
dimer to ferromagnetic phase. 

(ii) Generalized AKLT model, defined by F = 0, i.e. 
b = in the MP wave function (so that only triplets on 
the rungs occur): 

F = 0, x finite (17) 
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When choosing a convenient scale factor to render J e f t = 
1 this corresponds to the following Hamiltonian: 

J = 4/3 - 2x/3, Ji = J 2 = Ja = 5 /6, 

K = K' = 2/3, K" = Q, C = 3/4 + a;/3. (18) 

This is essentially the Hamiltonian of Ref . , but with- 
out requiring explicitly the coupling of two spins S = ^ 
into a triplet, a result that could be easily obtained di- 
rectly. The AKLT model in a strict sense is obtained for 
x — > +oo and a; FX 2 -) — » 0. 

(iii) For xF(x) = 1 another class of hamiltonians 
without biquadratic terms is obtained; the condition A2 > 
leads to the restriction x > §. A nontrivial case is 
x = i, which leads to the presence of only one diagonal 
exchange interaction, equivalent to a chain with NN and 
NNN exchange, 

C = 7/12, J = J = 2/3, J^J[ = J^ = -1, 

J 2 = K = K' = K" = 0. (19) 

This is a S = \ chain with alternating ferro- and anti- 
ferromagnetic exchange and ferromagnetic NNN interac- 
tions, a frustrated Heisenbcrg chain. As in (i), the eigen- 
value A2 vanishes and the energy of the fully polarized 
ferromagnetic state coincides with that of the MP ground 
state (||). Thus we have obtained two exact points on the 
line where the chain with NN and NNN interactions |L3| ] 
undergoes a first order quantum phase transition to the 
ferromagnetic state. 

(iv) A family of models with exact MP ground states 
connecting the Majumdar-Ghosh and the AKLT limits is 
obtained from e.g. 



F(x) = 1/(1 + xf. 



(20) 



Since the Majumdar-Ghosh and dimer chains have iden- 
tical ground states, Eq. (^o|) explicitly demonstrates the 
possibility to transform the Hamiltonian of the dimer 
S = ^ chain continuously to the Hamiltonian of the 
AKLT chain without any singularity (or quantum phase 
transition) in the ground state. 

This provides an exact proof that from the point of 
view of theoretical concepts no essential difference exists 
between the Haldane gap and the dimer gap and that the 
gapped S = j models listed in the introduction belong 
to the same phase as the Haldane S = 1 chain. 
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FIG. 1. Generalized spin-| ladder with all possi- 
ble isotropic exchange interactions between nearest-neighbor 
rungs. Ovals show spin pairs whose states are involved in the 
elementary matrix gt defined by (^). Couplings are described 
in the text, see (0). 
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